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A long and slender finger can serve as a simple "test bed" for different phase ordering models. 
In this work, the globally-conserved, interface-controlled dynamics of a long finger is investigated, 
analytically and numerically, in two dimensions. An important limit is considered when the finger 
dynamics are reducible to the area-preserving motion by curvature. A free boundary problem for 
the finger shape is formulated. An asymptotic perturbation theory is developed that uses the finger 
aspect ratio as a small parameter. The leading-order approximation is a modification of "the Mullins 
finger" (a well-known analytic solution) which width is allowed to slowly vary with time. This time 
dependence is described, in the leading order, by an exponential law with the characteristic time 
proportional to the (constant) finger area. The subleading terms of the asymptotic theory are also 
calculated. Finally, the finger dynamics is investigated numerically, employing the Ginzburg-Landau 
equation with a global conservation law. The theory is in a very good agreement with the numerical 
solution. 

PACS numbers: 64.75. +g, 05.70.Np, 05.70.Ln 
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I. INTRODUCTION 

This work is motivated by the recent developments 
in phase-ordering theory. Phase ordering is the growth 
of order from disorder via domain growth and coars- 
ening. As phase ordering systems are strongly non- 
linear and highly disordered, their theoretical descrip- 
tion remains challenging jlj. To get insight, it is of- 
ten useful to consider the coarsening dynamics of sim- 
ple test objects. One such object is a spherical droplet 
of the "minority" phase. It serves as the building 
block of simplified phase-ordering theories for locally con- 
served (bulk-diffusion-controllcd) (2| and globally con- 
served (interface-controlled) || systems. 

Being too a simple object, a spherical droplet may 
not give sufficient insight. For example, it's dynam- 
ics (shape-preserving shrinking or expansion) look the 
same in both locally, and globally conserved systems. 
On the other hand, there are important differences in 
the phase-ordering behavior of locally and globally con- 
served systems. In both cases, an ensemble of droplets 
exhibits Ostwald ripening, and the corresponding mean- 
field theories, due to Lifshitz and Slyozov || and Wag- 
ner respectively, predict dynamic scaling behavior of 
the droplet size distribution (with different dynamic ex- 
ponents). However, for finite volume fractions, the dif- 
ference in the type of conservation law leads to a dif- 
ferent role of correlations. Correlations between neigh- 
boring droplets are much more important in the locally- 
conserved systems [Q, than in the globally-conserved 
ones Even bigger differences have been found in 

the phase ordering dynamics of locally and globally con- 



served systems with long-range correlations. Coarsening 
of fractal clusters shows dynamic scale invariance and 
"normal" scaling in the case of global conservation ||, 
and breakdown of scale invariance and anomalous scal- 
ing in locally conserved systems . 

It is natural to attribute these differences to a basic 
difference in the character of transport: in globally con- 
served systems the transport is uninhibited by Laplacian 
screening effects, typical for locally conserved diffusion- 
controlled systems. To study this basic difference in a 
simple setting, the dynamics of a single long and slen- 
der finger is often considered. It has been found recently 
H that, in the locally-conserved, diffusion-controlled sys- 
tem, the finger acquires a dumbbell shape and shows 
non-trivial dynamic scalings for the finger length and the 
"ball" size, while the initial finger width remains almost 
constant until a late stage of the dynamics. Further- 
more, a finger-shaped domain (a "needle") served as a 
test object in still another conserved coarsening system: 
the one controlled by edge diffusion ||. Finger- like ob- 
jects appear naturally in two-dimensional simulations of 
dewetting ]ic[| , etc. In this work we investigate the glob- 
ally conserved dynamics of a long slender finger. An 
additional motivation for studying the finger dynamics 
is a recent observation that, at a late stage of coarsen- 
ing of fractal clusters, the cluster morphology shows long 
branches, or fingers, in both locally |?]||[LlJ , and globally 
HIH conserved systems. 

Here is an outline of the rest of the paper. In Section 2 
we briefly review the Ginzburg-Landau equation with a 
global conservation law (a phase field model for globally- 
conserved phase ordering) and its general sharp interface 
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formulation in two dimensions. Under a certain condi- 
tion (that will be elucidated) this formulation is reduced 
to the area-preserving dynamics by curvature. In Section 
3 we formulate a moving boundary problem for the fin- 
ger dynamics and develop a perturbation theory which 
enables us to obtain an intermediate asymptotic solu- 
tion for the finger shape. In Section 4 we return to the 
globally-conserved Ginzburg-Landau equation and inves- 
tigate the finger dynamics numerically. Section 5 includes 
a brief discussion. 



II. GLOBALLY CONSERVED PHASE 
ORDERING: A PHASE FIELD MODEL AND ITS 
SHARP-INTERFACE LIMIT 

Globally-conserved phase ordering can be described by 
a simple phase field model Jl5|-|l5|] . This model represents 
the dynamics in terms of a simple gradient descent 



du _ 5F 
dt Su 

with the Ginzburg-Landau free energy functional 



F[u] 



[(1/2)(VV) 2 + V{u) + Hu] dr 



(1) 



(2) 



and a double- well potential V(u) = (1/4)(1 - u 2 ) 2 . The 
effective "magnetic field" H = H(t) varies in time so as 
to impose global conservation law 



< u(r, t) >= 



: f u(r, t) dr — const . 

j xL y J 



(3) 



Here u(r, t) is the order parameter field, L x 3> 1 and 
L y 1 are the linear dimensions of the system (a two- 
dimensional rectangular box), and the integration is car- 
ried out over the entire box. Eqs. ([[])- (Q) yield the non- 
local Ginzburg-Landau equation 



du 
0* 



= V 2 u + u-u 6 



< u ■ 



u 3 > 



(4) 



Either no-flux, or periodic boundary conditions are as- 
sumed. 

For a wide class of initial conditions this coarsening 
system segregates, at late times, into large domains of 
"phase 1" and "phase 2" separated by thin domain walls 
(which width is of order unity) ||[l3]|l(|. Correspond- 
ingly, a sharp interface theory can be developed for these 
late times H. At these times the magnetic field H(t) is 
already small, H(t) <C 1, and slowly varies with time. 
The phase field in the phases 1 and 2 is almost uniform 
and rapidly adjusts to the current value of H(t), so that 
u = -1 - H(t)/2 and 1 - H{t)/2, respectively. The nor- 
mal velocity of the interface is given by M 



v n {s,t) = k{s, t) - ■^ H ( t ) ■ 



(5) 



where k is the local curvature and s is the coordinate 
along the interface. The positive sign of v n corresponds 
to the interface moving towards phase 1, while k is posi- 
tive when the interface is convex toward phase 2. 
The dynamics of H(t) is determined by || 



H(t) 



4A(t) 

L X Ly 



k(s, t) 



V2 



H(t) 



where 



k(s, t) 



A(t) 



k(s, t) ds , 



(6) 



(J) 



and A(t) is the total interface length. Equations (||) and 
(^) make a complete set and provide a general sharp- 
interface formulation to this problem. In some cases this 
formulation can be further simplified |5|,[l3| . Let us com- 
pute the rate of change of the domain area, A: 



A(t) = <p v n (s,t)ds = A(i) 



K{s,t)--j=H(t) 



(8) 



If the two terms in the right hand side of Eq. (||) balance 
each other, 



V2- 



H(t)~—K(s,t), 



(9) 



the domain area remains approximately constant. Then, 
using Eqs. (pi) and (|9|), we obtain 



v n (s,t) = n(s,t) - n(s,t) . 



(10) 



Model (10) is known as area- preserving flow by curva- 
ture p7|. Dynamics ( |l0| ) reduce the interface length of 
the system |L3 17 1, leaving the areas of each of the phases 
constant. The only stable steady state of Eq. (jll]) (not 
attached to the boundaries of the system) is a single 
perfectly circular domain p3| . Under what condition is 
model ( [l0| ) a good approximation to the more general 
sharp-interface model (|5|) and (0)? Consider Eq. (|^) 
which includes the same combination k — (3/V2)H(t). 
Let us treat the term in the right hand side of Eq. (|^) 
perturbatively, putting H(t) — H^°'(t) + h(t), where 
H(°\t) = (y/2/3)K is the leading term, and h(t) is a 
subleading term. Keeping terms up to order h(t) we ob- 
tain: 



h(t) 



y/2L x Ly 
12A 



ffW(t) 



L X Ly 

18A 



(11) 



If there is only one (single-connected) domain, k — 
— 27r/A. Therefore, 



h(t) 



irL x Ly A 
9 A3 



(12) 



Requiring \h{t)\ < |F (0) (i)|, 

we obtain the "area- 
preservation" criterion 



2 



L x L y \A\ 
A 2 



«1. 



(13) 



In Section 3 we will investigate the area-preserving dy- 
namics of a long slender finger, and verify criterion Jl3| ) 
a posteriori. 



III. AREA-PRESERVING FINGER DYNAMICS: 
A PERTURBATION THEORY 

Let the initial condition for dynamics ( |To| ) represent 
a long and narrow rectangular bar with length 2a and 
width 2A . We will see that, with time, the bar evolves 
into a finger-shaped object. We place the origin of a 
Cartesian coordinate system in the finger's center and 
denote by x the coordinate along the finger and by y — 
y(x, t) the instantaneous location of the finger boundary. 
Because of the symmetry with respect to the ic-axis, we 
will be interested only in the upper boundary dynamics: 
y(x, t) > 0. In terms of y(x, t) the local curvature is given 
by @: 



k(s, t) 



JJx 



(l+^) 3/2 



where y x = dy(x,t)/dx. The average curvature n(s,t) is 

— r — 27T — 7T 

K(S, t) 



A 



2Jo° (t) (l + y 2 x )^ 2 dx 



where a(t) is the time-dependent position of the finger 
tip. Since 



^ = (1 

dt 1 



2U/2, 



dynamics ([L0|) imply the following evolution equation for 
y(x,t): 



dy_ 

dt 



t(1 + J^) 



2U/2 



1 



Vl 2j a{t \l + yl)^dx 



(14) 



a(t) is defined by the boundary conditions y(x = a, t) = 
and y x (x — a.t) = — oo. An additional boundary condi- 
tion follows from the symmetry with respect to the y-axis: 
y x (x = 0, i) = 0. Equation (14), a nonlinear integro- 
differential equation in partial derivatives, describes a 
free-boundary problem. It is easy to check that the dy- 
namics ( |T4| ) preserve the finger area A. 

It can be proved Jl^ j that any convex domain evolving 
by Eq. ( p~0| ) [and therefore by Eq. (|i"4|)] finally becomes 
a circle with the same area A. We are interested in the 
finger dynamics at times much smaller than those needed 
for approaching this equilibrium (but still large enough 
so that details of the initial shape are forgotten). The 
corresponding strong double inequality will be presented 
at the end of this Section. 



As long as the ratio of the finger width to its length 
is small, we can use it as the small parameter of our 
theory. Indeed, near the finger tip the first (local) term 
in the right hand side of Eq. ( [bi] ) is of the order of the 
inverse finger width, while the second (nonlocal) term 
is of the order of the inverse finger length. Therefore, 
the nonlocal term can be treated perturbatively. The 
unperturbed equation, 



dy 

dt 



y.i 



(15) 



was obtained more than 40 years ago by Mullins in the 
context of motion of grain boundaries 19 1 . It corresponds 
to a non-conserved motion by curvature: v n (s,t) = k. 
We notice immediately that the perturbed problem, Eq. 
(|l4|), possesses an integral of motion (the finger area) 
which is absent in the corresponding unperturbed prob- 
lem. This situation is unusual: more often than not per- 
turbations destroy integrals of motion pertinent to the 
unperturbed system. This unusual property will be ex- 
ploited in the following, as we will require area conserva- 
tion in each order of the perturbation expansion. 

Mullins obtained a one-parameter family of 

traveling-wave solutions of Eq. ( |l5| ) describing a half- 
infinite constant-width finger retreating by its tip's cur- 
vature: 



y = Y(x, t) = - arccos [exp (— c£)l , 
c 



(16) 



where £ = a(t) — x, a(t) = — ct is the time-dependent 
position of the finger tip, and the constant speed of re- 
treat c is the parameter of this family of solutions. The 
finger's half- width A = 7r/(2 c) is constant. We will show 
that this solution (with two important modifications) can 
serve as the leading (or zero) order approximation of 
our perturbation theory for the area-preserving finger dy- 
namics. First, we introduce a (slow) time-dependence in 
the parameter c. Second, we "stick together" two iden- 
tical (very long but finite) half-fingers. The resulting 
ansatz is the following: 

Y(x, t) = -j- arccos [exp (-c(t) £)] , (17) 
c(t) 

where £ = a(t) — \x\, and a(t) > is yet unknown. The 
finger half- width A(t) = n/2c(t) is now time-dependent; 
its growth with time makes up, in the area conservation, 
for the finger shortening. As a function of x, ansatz ( |r^ ) 
is continuous everywhere. Its ^-derivative is not contin- 
uous at x = 0. However, as far as e(i) = A(t)/a(t) is 
very small: e(t) <C 1, the x-derivatives at x = from left 
and right are of order of exp(— 1/e), that is exponentially 
small. This and other exponentially small effects will be 
neglected throughout the paper. Because of the symme- 
try with respect to the y-axis, we will consider only the 
region < x < ait). The perturbation expansion for 
y(x,t) is 



3 



y(x,t) =Y(x,t)+Sy(x,t) + 



(18) 



where 



where 5y(x, t) is the subleading (or first order) term. In- 
serting Eq. ( |l7| ) into Eq. ( |l4| ) and keeping only the zero 
order terms, we obtain 



= -c= -7r/(2 A). 



(19) 



The evolution equation for A(t) follows from the require- 
ment of area conservation in the zero order. Calculating 
the area under the Y{x, t) profile, we obtain: 



»(*> . . 2A 2 ln2 
Y(x, t) dx = aA 



0(exp(-l/e)). 



(20) 



One can see that the finger area in the zero order is simply 
A = A a A. Demanding A — and using Eq. @, we 
obtain: 



A = 7r/(2 a) 



(21) 



We will see later that the same equation follows from the 
analysis of the first order correction. The solution of the 
zero order equations (E9h and (Ell) is 



and 



a (0) (i) = a exp (-7rf/2a A ) , 



A^°\t) = A exp (7rf/2a A ) 



(22) 



(23) 



Therefore, in the zero order, the finger shortens (and its 
width grows) exponentially with time. The characteristic 
growth time is of order of the finger area. The zero order 
solution yields a criterion for the validity of the pertur- 
bation theory. The aspect ratio of the evolving finger 
should be very small which leads to 



t <C a A In 



Of) 



(24) 



This condition sets an upper limit for the times for which 
our perturbation expansion is valid. 

We now turn to calculating the small correction 
5y(x, t). Since the nonlocal term of Eq. (Jl4|) is already of 
order e, we evaluate it on the zero-order solution Y(x, t) 
which yields 7r (2a) _1 (l + Y 2 ) 1 / 2 . Introduce a new vari- 
able u: 



u(x,t) = exp(-c(i)£) . 



(25) 



Substituting Eq. ( |l8| ) into Eq. ( |l4"| ) and linearizing, we 
get a linear partial differential equation for 5y(x,t): 



Sy t = (1 - u 2 ) Sy xa: 
(c + d) u 



2cu Sy x 
c 



(l-u 2 ) 1 / 2 & 



q(u) 



2a(l-w 2 ) 1 /2 



(26) 



q(u) 



u In u 



{l-u 2 ) 1 / 2 



arccos u . 



(27) 



It is convenient to go over to new independent vari- 
ables £ and t. We define 5y(x,t) = /(£, i) and obtain 
Syt = ft + fl/f- As the rj-dependence in the new vari- 
ables is slow, the term f t can be neglected in this order 
of the perturbation scheme. Therefore, we are left with 
an ordinary differential equation (the slow time enters as 
a parameter): 



(1 - u 2 )fe + c(2u 2 + l)f, = , 1 (C+ .t ) r /9 - 4f<?('" > 



(l-u 2 ) 1 / 2 c 2 ' 

7T 

~2a{l-u 2 ) 1 / 2 ' 



(28) 



Going over from £ to u and defining /(£) = g(u) we 
obtain 



(1 - u 2 )g uu - 3ug u 



(c + a) 



1 



c q(u) 



c 2 uil-u 2 ) 1 / 2 c 4 u 2 

7T 



2ac 2 u 2 {\~u 2 ) 1 / 2 



(29) 



Define <I>(ii) = g u (u). The general solution of the homo- 
geneous equation for <E>(it) is CQq(u), where 



$oW - (1 - u 2 )- 3 / 2 , 



(30) 



and C is an arbitrary constant. Therefore, we look for 
the general solution of Eq. ( p9| ) in the form <&(w) = 
C(u)Qq(u), where C(u) is an unknown function. Substi- 
tuting this into Eq. ( p9| ) and integrating, we obtain 

C(u) = £±£ Inu- A a(u) + + a {t) , (31) 



where 



2ac 2 i 



s(u) — — In 2 u H — arccos 2 u — In u 
K ' 2 2 



{l-u 2 ) 1 ' 2 



arccos u 



(32) 



and a(t) is an arbitrary function of time. Going back to 
g(u), we can write the first order correction to the finger 
shape as 



9{u) 



C{u)§ (u) du + /3(t), 



(33) 



where (3{t) is another arbitrary function of time. The 
boundary conditions for g(u) are g(0) = and g(l) = 
0. The first boundary condition gives P(t) — 0, while 
the second boundary condition should be used to find 
an equation for d in the first order. Before doing this, 
an analysis of the integrals in Eq. pS] ) should be made. 
To prevent divergence of g{u) at u — 1 we must require 



4 



a(t) = — 7r/(2ca). Similarly, to prevent divergence of 
g(u) at u — we must require c/c 2 = —1/a. The latter 
equation is equivalent to Eq. obtained in the leading 
order from the requirement of area conservation. 

The equation for d in the first order of the perturbation 
theory is obtained from the boundary condition g(l) = 
which requires calculating the integral in Eq. (|3^) from 
u = to u = 1. This calculation yields 



3 In 2 



7T 

2A 



3 In 2 



(34) 



It is now possible to obtain an explicit expression for 
g(u). An integration from to u in Eq. (33) gives: 



9{u) 



1 



w\(u) — arcsin uln 



2 n l e[l + (l-w 2 ) 1 /2 




(35) 



where 



and 



2 u In it In (8u 1 / 2 /e) + it arccos 2 u + n(l — u) 



2(l- U 2)l/2 



u> 2 (m) = 2 



arcsm u 



du , 



(36) 



(37) 



In the absence of a conservation law, we would have to 
go to the second order of the perturbation scheme to 
obtain the equation for c in the first order. In our area- 
preserving perturbation scheme this equation can be ob- 
tained already in the first order. Notice that the first- 
order correction obeys the following conservation law: 



£ / ° 6y(x,t)dx = a[f(t = a,t)-f(Z = 0,t)} = 0, 



(38) 



d 
dt 



where the last equality results from the boundary condi- 
tions. This relation combined with the requirement that 
the finger area is preserved leads to 



d_ 

dt 



a(t) 



Y(x,t)dx = 0. 



(39) 



It should be stressed that functions a(t) and A(t) enter- 
ing Eq. (|39|) now include first order corrections. Now, 
using Eqs. ((2^), (^ ) and (^9| ) and keeping terms only 
up to the first order in e, we obtain 



7T 

2a 



A In 2 



(40) 



Equations ([34|) an d ( j40| ) are the equations for d and 
for A in the first order of the perturbation scheme. This 



set of equations is autonomous and therefore integrable. 
In the framework of the perturbation theory it is more 
consistent to solve the equations perturbatively by lin- 
earization. To do this, we exploit the conservation law 
( p9| ) which, together with Eq. (^0|), allows us to express 
A(t) as a function of a(t) or vice versa. The expression 
obtained is then linearized with respect to A/a and in- 
serted into Eq. (||) (or @). This yields the following 
equations: 



4 In 2 



i In 2 



2a A 



A = 



7rA 

2a A 



A In 2 A In 2 A 3 In 2 



a 2 A 2 



(41) 



(42) 



We are looking for solutions in the form of a(t) = (t) + 
a^(t) and A(t) = A^{t) + A«(t), where a^(t) and 
A^(t) are first order corrections. Solving the resulting 
linear equations, we obtain: 



* (1) (*) = 



8A In 2 



sinh 



(- 



nt 



n 

tin 2 
a 



■ exp 



\2aoA 

7rf 



2a A 



(43) 



A«(t)= 4A « ln2 



irao 
A tln2 



• exp 



exp 



nt 
aoA 



sinh 



Tit 



V2a A 



V 2a A 



(44) 



We checked that when solving Eqs.(J34|) and (^) numer- 
ically, the results differ from the approximate analytical 
results (^3|) and (^4|) by less than 1%. Finally, requiring 
a<°>(t) > a {1 \t) and A^(t) > A«(t) and using @, 
(p3|), (p3|), and (Q), one can see that the linearization 
procedure is valid for times obeying the same strong in- 
equality (^4|). Therefore, the linearization procedure is 
consistent with the perturbation theory. 

The validity of our perturbative solution is limited by 
not too long times, see inequality (p4]). On the other 
hand, one should wait until irrelevant details of the ini- 
tial bar shape are forgotten, and the solution approaches 
Y{x,t). Typically, it takes time needed for the bar tip 
to pass the distance of order of the characteristic tip size 
1/c ~ A. The tip speed is of order 1/A, so the "waiting 
time" of the theory can be estimated as Aq. Our per- 
turbative solution represents, therefore, an intermediate 
asymptotics, valid for 



Al < t < a A In ^~ 
Aq 



(45) 



Finally, we should check the area-preservation criterion 
(|l3|). In the leading order we have: A(t) ~ a^°^(t) and 



5 



A(t) ~ dW(t) - 1/A(°)(t). Using (H) and © in flij) 
we obtain: 



i < a A ln(/a ) , 



(46) 



where / = (a A ) / '(L x L y ) is the area fraction occupied 
by the finger. 



IV. NUMERICAL SOLUTION 

In order to test the predictions of the perturbation the- 
ory we performed numerical simulations with the nonlo- 
cal Ginzburg-Landau equation (Q) (with no-flux bound- 
ary conditions), taking a long and narrow rectangular 
bar as the initial condition. We used an explicit Euler 
integration scheme to advance the solution in time, and 
second order central differences to discretize the Laplace 
operator. A mesh size of Ax = Ay = 0.5 was found 
sufficient for an accurate resolution of the interface. A 
time step of At = 0.05 was required for numerical sta- 
bility. The bar width 2Ao = 50 was chosen to guarantee 
area conservation with a good accuracy [see Eq. (ftq)]. 
The simulation was carried out up to a time tf = 9 • 10 4 
which was long enough to distinguish between an expo- 
nential and linear dynamic behavior for a(t) and A(t). 
The bar length 2ao = 8 • 10 3 was chosen, so that the 
parameter e(t) — A(t)/a(t) was sufficiently small even 
at the end of simulation (when it reached 0.1). The bar 
was placed in the center of a rectangular box 8020 x 220. 
Because of symmetry with respect to the x and y axes 
only the quadrant x > , y > was actually simulated. 
We checked that the finger area was conserved with an 
accuracy better than 1%. 

The time dependence of the finger half-length a(t) ob- 
tained from the numerical simulations is shown in Fig. 
1. The same graph also shows the theoretical values of 
a(t) in the zeroth order Wf\t), Eq. (|22|)1, and first order 
[d°\t) + aW(f), Eq. (||)] of the perturbation theory. 
A good agreement between the numerical and analyti- 
cal results is obtained already in the zeroth order. This 
agreement is improved further by the first order correc- 
tion. For example, at the final time of the simulation 
tf =9 - 10 4 , there is a 9% difference between the numer- 
ical value of a(t) and a^°'(t), and only a 1.4% difference 
between the numerical value of a(t) and (t) +a^\t). 

The numerically found half- width of the finger A(t) 
is shown in Fig. 2. Its comparison with the zero-order 
values A(°)(t) [Eq. (|23|)] and with the first-order values 
A (a \t) + AW(t) [Eq. (j4§)] is also shown. Again, there 
is a good agreement already in the zeroth order, and this 
agreement is further improved by the first order correc- 
tion. At the final time of the simulation, the numerical 
value of A(t) differs by 6% from A^ ^ (t), and by less than 
2% from A (0) (i) + A^(t). 
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FIG. 1. A semi-log plot of the finger half-length a versus 
time. The solid line represents a{t) found numerically. The 
dashed line shows the zero-order analytic result a^°'(i). The 
dotted line is the first-order analytic result a^°\t) + a'^{t). 




2x1 4 4x1 4 6x1 4 8x1 4 
t 

FIG. 2. A semi-log plot of the finger half-width A versus 
time. The solid line is A(t) found numerically. The dashed 
line shows the zero-order analytic result A^°'(i). The dotted 
line is the first-order analytic result A (0) (i) + A (1) (i). 

Fig. 3 shows the finger shape y(x, t) found numerically, 
at the final time of the simulation tt = 9 • 10 4 (notice 
the different scales in the x and y axes). Also shown 
are the following zero- and first order theoretical results 
evaluated at t = tf: 
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FIG. 3. The finger shape y(x,tf) (solid line). The dashed 
line is the function Y m (x,t f ) defined by Eq. (ff7|), the 
dashed-dotted line is the function (fji|), and the dotted line 
is the function (^). All the functions are evaluated at 
t = tf = 9 • 10 4 . The dotted line is indistinguishable from 
the solid line. 

Y^(x,t) = -^arccos [exp (- C (°>(i)£ (0) )], (47) 
where = a^(t) - x, and 
Vs{x, t) = y— ■ arccos [exp (-Cs (*)£,)] + Sy(x, t) , (48) 

C s (t) 

where £ s = a s (i)-x, a s = a^+a^\ and c s = c^+cPK 
In addition, shown is the function Y num (x,t) defined by: 

Y num (x, t) = 5—— arccos [exp {-c num (t)^ num )] , (49) 

Cnum \t) 

where £ num = "numW and a, n um{t) and c num (t) are 
the numerical values of a(t) and c(t), respectively. The 
function Y num (x, t) is also evaluated at t = tf. The com- 
parison with the zero-order prediction (x,tf) shows 
a good agreement even at this relatively late time, when 
A/a ~ 0.1. The agreement is further improved by 
the first-order theory. The comparison with Y num (x,tf) 
shows an excellent agreement, implying that the sublead- 
ing term Sy(x,t) is very small compared to the leading 
term Y(x,t). Fig. 4 which shows the theoretical result 
for 5y(x,t) at t = tf, calculated by using Eqs. (|35|)- 
(|37|), confirms the later conclusion. We should add that 
a very good agreement was found between the theoreti- 
cal Sy(x,t) and the numerical 5y(x,t), which is obtained 
by subtracting Y num (x,t) from the numerical result for 
y(x,t). 



200 400 600 800 1000 
X 

FIG. 4. The subleading term Sy(x,tf) at the final time of 
the simulation, as obtained by the perturbation theory [Eqs. 
©-©]• 

V. CONCLUSIONS 

We investigated the dynamics of a long slender finger- 
shaped domain undergoing globally-conserved (interface- 
controlled) coarsening. We worked in the parameter 
region where these dynamics can be reduced to area- 
preserving motion by curvature. The nonlinear moving 
boundary problem for the finger shape was solved by us- 
ing an asymptotic perturbation theory which employed, 
as the zero-order solution, a straightforward modifica- 
tion of the Mullins solution for a finger retreating by 
the curvature of its tip. Both the leading, and the sub- 
leading terms of the solution were calculated analytically 
and verified by a numerical solution of the correspond- 
ing phase-field equation: the Ginzburg-Landau equation 
with a global conservation law. 

The coarsening dynamics of the finger look quite dif- 
ferently from those observed in the case of a locally con- 
served (diffusion-controlled) system ||. The finger pre- 
serves its simple shape, while its characteristic length 
and width are changing on the same time scale. There- 
fore, the global character of transport (uninhibited by 
Laplacian screening effects characteristic of the locally 
conserved system) manifests itself already in the simple 
setting of finger dynamics. 

ACKNOWLEDGEMENTS 

This work was supported in part by the Israel Science 
Foundation, administered by the Israel Academy of Sci- 
ences and Humanities. 



7 



A. J. Bray, Adv. Phys. 43, 357 (1994). 

I.M. Lifshitz and V.V. Slyozov, J. Phys. Chem. Solids 
19, 35 (1961). 

C. Wagner, Z. Electrochem. 65, 581 (1961). 
M. Marder, Phys. Rev. A 36, 858 (1987). 

B. Meerson and P.V. Sasorov, Phys. Rev. E 53, 3491 
(1996). 



A. Peleg, M. Conti, and B. Meerson, :ond-mat/000722C 
M. Conti, B. Meerson, and P.V. Sasorov, Phys. Rev. Lett. 
80, 4693 (1998). 



M. Conti, B. Meerson, and P.V. Sasorov, e-print cond 



mat/9912426 



R. Thouy, N. Olivi-Tran, and R. Jullien, Phys. Rev. B 
56, 5321 (1997). 

J. Koplik and J.R. Banavar, Phys. Rev. Lett. 84, 4401 
(2000). 

T. Irisawa, M. Uwaha and Y. Saito, Europhys. Lett. 30, 
139 (1995). 

T. Irisawa, M. Uwaha and Y. Saito, Fractals 4, 251 
(1996). 

J. Rubinstein and P. Sternberg, IMA J. Appl. Math. 48, 
249 (1992). 

C. Sire and S.N. Majumdar, Phys. Rev. E 52, 244 (1995). 
A.D. Rutenberg, Phys. Rev. E 54, 972 (1996). 
A.S. Mikhailov, Foundations of Synergetics. 1. Dis- 
tributed Active Systems (Springer- Verlag, Berlin, 1993). 
M. Gage, Contemp. Math. 51, 51 (1986). 
R. Courant, Differential and Integral Calculus, Vol. /(In- 
terscience Publishers, New York, 1956). 
[19] W.W. Mullins, J. Appl. Phys. 27, 900 (1956). 



8 



